ITP-UU-12-22 
SPIN-12-20 
TCD-MATH-12-05 
HMI-12-02 

Prepared for submission to JHEP 



Exceptional Operators in J\f = 4 super Yang-Mills 



Gleb Arutyunov,"' 1 Sergey Frolov 6,1 and Alessandro Sfondrini a 

a Institute for Theoretical Physics and Spinoza Institute, Utrecht University, 

Leuvenlaan 4, 3584 CE Utrecht, The Netherlands 
h Hamilton Mathematics Institute and School of Mathematics, 

Trinity College, Dublin 2, Ireland 

E-mail: G.E.Arutyunov@uu.nl, f rolovsOmaths . ted. ie , A.Sfondrini@uu.nl 

Abstract: We consider one particularly interesting class of composite gauge-invariant op- 
erators in M = 4 super Yang-Mills theory. An exceptional feature of these operators is that 
in the Thermodynamic Bethe Ansatz approach the one-loop rapidities of the constituent 
magnons are shown to be exact in the 't Hooft coupling constant. This is used to propose 
the mirror TBA description for these operators. The proposal is shown to pass several 
non-trivial checks. 



Correspondent fellow at Steklov Mathematical Institute, Moscow. 



Contents 



1 Introduction and summary 1 

2 Bethe-Yang equations and singular rapidities 6 

2.1 Singular rapidities in the one- loop Bethe Ansatz 7 

2.2 All- loop Bethe-Yang equations and their breakdown 8 

3 Twisted TBA 10 

3.1 Analytic structure of Y-functions 10 

3.2 Wrapping corrections for L = 6 at 0{g 12 ) 11 

3.3 Comments on the g 14 correction 12 

4 TBA with exceptional rapidities 14 

4.1 TBA equations 15 

4.2 Exact Bethe equations 17 

4.3 Leading TBA correction up to g w 20 

4.4 Next-to-leading TBA correction at g 12 21 

5 Conclusions 24 

6 Appendices 26 

6.1 Twisted transfer matrices and relating twist to a 7-deformation 26 

6.2 Twisted Y-functions and their analytic properties 27 

6.3 TBA equations for the twisted theory 28 

6.4 Linearized TBA and exact Bethe equations for L = 6 31 



1 Introduction and summary 

The aim of this work is to provide the mirror TBA description of one particularly interesting 
class of composite gauge-invariant operators in planar N = 4 super Yang-Mills (SYM) 
theory and thus to further advance understanding of the planar AdS/CFT [1] spectral 
problem. 

The operators we are interested in belong to the so-called su(2) sector of the M = 4 
SYM and they are eigenstates of the one-loop dilatation operator having the following 
explicit form [2] 

L-4 

O l = J2(-iy tx(xX Z { X Z L - { -^ . (1.1) 

i=l 

Here X and Z are complex scalars of N = 4 SYM and L > 6 is an even number. 



- 1 - 



Our special interest in this class of operators is motivated by the following. At one 
loop operators from the su(2) sector can be identified with excitations of the XXX Heisen- 
berg spin chain [3]. From this point of view, the operators above represent three-particle 
(magnon) states, and the simplest of them is an excitation of the spin chain of length L = 6. 
Diagonalizing the Heisenberg Hamiltonian for this case, one finds the corresponding eigen- 
value to be where A is the 't Hooft coupling. Thus this state is in the spectrum of the 
XXX model and the same conclusion holds for all Ol- However, trying to describe these 
states by solving the corresponding Bethe Ansatz equations one encounters a problem - 
the magnons must have their rapidities Uj at distinguished positions in the complex plane, 
namely at 1 — i, 0,i [4, 5]. As a result, the scattering matrices entering the Bethe Ansatz 
are singular and the energies of such states are ill-defined 2 . This problem is, of course, 
well known and one natural way to cure it is to introduce a regularization by means of a 
twist, which we call (ft. In the gauge theory twisting can be linked to the Leigh-Strassler 
deformation of M = 4 super Yang-Mills theory [6] dual to strings in the Lunin-Maldacena 
background [7] with a real deformation parameter and their nonsupersymmetric general- 
izations [8]. In this physical theory the limit (ft — > can be taken without any problem. In 
the Bethe Ansatz approach one first computes the energy of Ol for finite (ft and then takes 
finding the same result as from the direct diagonalization of the Hamiltonian. 

Also, having rapidities of two magnons at singular points ±i can be related to the fact 
that Ol is a mixture of operators where two fields X are stuck together. In the terminology 
of [9] two magnons form an infinitely tight bound state. We will have to say more about 
the nature of this bound state later. 

Obviously, at one loop introduction of a twist is just a minor feature which distinguishes 
Ol from other operators. Going to higher loops reveals more dramatic differences. To 
analyze the states corresponding to Ol at higher loops, we can try to employ the all-loop 
asymptotic Bethe Ansatz [2], which is also referred to as the Bethe- Yang equations. In 
addition to the twist (ft the Bethe- Yang equations depend on the coupling constant g which 
we identify with the effective string tension related to A as g = ^ . Expanding the Bethe- 
Yang equations in powers of g and starting from the one-loop rapidities 0, ±i, one can find 
a formal power series solution for Uj with coefficients depending on (ft. As expected, nothing 
special happens until one reaches the first wrapping order. However, at the first wrapping 
order, g 2L , one discovers that the limit (ft — > is singular and the corresponding energy 
diverges as (ft approaches zero. This behavior should be contrasted to that of regular states 
{e.g. Konishi): the latter do not even require the introduction of a twist. On the other 
hand, from the point of view of the gauge theory we should do not expect any problem 
with taking 0—7-0 for operators of the type Ol- 

Certainly, the Bethe Ansatz is only asymptotic, that is it provides a correct description 
of the spectrum only up to the first wrapping order; the perturbative behavior of Ol serves 

1 Here Bethe roots are rescaled by a factor 1/2 in comparison to the XXX standard normalization. 
2 At one loop one can use Baxter's Q-operator to describe the corresponding states in terms of dual roots 
which lead to the well-defined energy. 
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as a clear confirmation of this fact. Hence, as for regular operators, we should expect that 
the TBA must give an adequate solution. 

We recall that the TBA approach, originally developed for relativistic theories [10], 
enables a computation of the ground state energy of a two-dimensional integrable model 
in a finite volume by evaluating the partition function of the accompanying mirror model 
[11]. In recent years the mirror TBA - a tool to determine energies of string states on 
AdSs x S 5 and correspondingly scaling dimensions of gauge theory operators - has been 
largely advanced [12]-[23] and generalized to include excited states [24]-[27]. Results derived 
from the corresponding TBA equations [28]- [31] show an agreement with various string 
[32]- [35] and gauge theory [36]- [38] computations, and also with Liischer's perturbative 
treatment [39]- [42]. 

Apparently, constructing the TBA equations for the states corresponding to Ol we 
might follow the same procedure as for regular states. This amounts to first building up an 
asymptotic solution with a finite twist 3 , analyzing its analytic properties and then using 
them to engineer the TBA equations [21]. However, the TBA equations constructed in 
such a way rely on the asymptotic solution which is valid only for A < which makes 
obscure how to take the limit <j> — > with A fixed. More precisely, for fixed <j> there always 
exists a critical value A cr = A cr (0) such that the Bethe-Yang equations have a well-defined 
solution for A < A cr (0) and no solution for A > A cr ((/>). Nevertheless, in perturbative 
treatment of the TBA this problem of order of limits can be overcome by considering first 
the expansion in powers of A and then taking the limit <f> ~ > in each term of the expansion. 
In this work we consider in detail the corresponding twisted TBA equations for Ol with 
L = 6. In fact, introduction of the twist results in the analytic behavior of rapidities and 
Y-functions very similar to that considered in [27], in particular, the complex rapidities ui 
and U3 of the second and third particle respectively, lie outside the analyticity strip, which 
is in between two lines running parallel to the real axis at — ~ and ~. Not surprisingly, 
the TBA equations for the state corresponding to Oq essentially coincide with that of [27]. 
By expanding these TBA equations up to A 6 , we then show that the TBA correction to 
the Bethe-Yang equations cancels precisely the divergent part of the asymptotic energy 
rendering therefore the limit <j> — > well-defined. For the energy 

£(6) 

at six loops (the first 

wrapping order) we then find 

, = /_ 84753 243 189 _ 567 \ 

V 1024 128 sv ; 64 sw 128 s y J yv> ) y 

We also provide a mechanism for a similar cancellation at higher orders of A. This in 
principle solves the problem of describing singular states in perturbative theory. It is quite 
remarkable that in spite of the fact that the TBA corrections make the energy of a state 
finite in the limit <p — > 0, the perturbative rapidities remain divergent in this limit. 

A veritable question is however how to describe singular states for finite A and what 
are the corresponding TBA equations. To answer this question, we again consider a state 
which contains only our three distinguished magnons. For large L such a state can be 



Introduction of a twist in the mirror TBA has been considered in the recent work [43]- [46] 
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viewed as a scattering state of a fundamental particle and a two-particle bound state with 
momenta ±ir. We put forward a conjecture that the one-loop rapidities are in fact exact for 
any value of A, and we use this conjecture to propose TBA equations for these states. In 
what follows we refer to these rapidities as exceptional. As a very non-trivial consistency 
check, we show that our conjectured TBA equations lead to the constraints 4 Y\ t (uj) = — 1 
which for regular states would have to be imposed as momentum quantization conditions. 
We compute the energy of the shortest operator of this type (of length L = 6) up to A 6 
and show that it perfectly agrees with the result obtained from the twisted TBA equations. 
We believe that the equality of energies computed from the TBA based on twisted and 
exceptional rapidities must hold to all orders in perturbation theory. 

Amazingly, in the approach based on the exceptional rapidities, the TBA corrections 
begin to contribute to the energy already at A 3 , and for a generic singular state at A i//2 , i.e. 
at half -wrapping. This behavior is consistent with the analysis of [11] where a two-particle 
bound state with the total momentum p larger than the critical value p CT has been studied. 
Indeed, the leading exponential correction to the energy of the bound state found from 
the Bethe-Yang equations is e~ qJ and the leading TBA correction is expected to be of the 
same order. Here q is used to parametrize the complex particle momenta Pi = \ J riq and 
P3 = | — iq with Re q > 0. At weak coupling p cr ~ tt — 2g, i.e. the momentum p = tt of the 
two-particle bound state we are interested in here exceeds the critical value. According to 
[11], in the limit g — > one has q = — log | + . . ., i.e. the leading TBA correction must be 
of the order e~ qJ ~ g J where J ~ L is large. 

The family of three-particle states corresponding to Ol is probably the only example 
of states which rapidities are known as exact functions of g. For this reason we call the 
operators Ol exceptional. In a sense these states are similar to the vacuum state for which 
one does not have the exact Bethe equations. Of course, the TBA equations for Ol are 
non-trivial and they are ultimately responsible for the non-trivial dependence of energy on 
the coupling constant. It would be very interesting to see whether Ol exhibit exceptional 
features also from purely field-theoretic point of view. 

In fact one can consider more general operators which include the three exceptional 
magnons as a building block [5]. In contrast to the exact rapidities of exceptional magnons, 
extra rapidities of such an operator are not rigid and have non-trivial A-dependence. The 
results of this paper allow one to readily construct the corresponding TBA equations. In 
a sense all such states can be viewed as a new sector of N = 4 SYM with exceptional 
operators playing the role of non-BPS vacuum states. 

Having established two TBA approaches to exceptional operators - the twisted one and 
the one based on the exceptional rapidities (both producing the same perturbative energies) 
- one can naturally wonder what is the relation between them. Apparently, they look rather 
different, in particular, in the twisted approach the perturbative rapidities are divergent in 
the limit <p — > 0. To clarify this issue, one can fix a value of A and look for the evolution 
of the rapidities Uj(4>) when the twist decreases from some finite value to zero. Inverting 
the function A cr ((/>), one finds a critical value of the twist (j) cv = cr (A). For (p > (/> cr (A) 

4 Here Yi„ denotes analytic continuation of the main Y-function Yj to the string region. 
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Figure 1. The picture of the z-torus at g = 0.1. Brown and green curves are the boundaries of 
the string and (anti-)mirror regions. They intersect at four points which correspond to the branch 
points on the string u-plane. Exceptional rapidities are u± = 0, 1*2,3 = Ti/g- Twisted rapidities 
are indicated by Vi and Wi. Rapidities 1*2.3 are located just a bit outside the analyticity strip as 
happens for e.g. L = 6, while u>2,3 are inside as for e.g. L = 8. Arrows indicate the conjectured 
motion of the twisted rapidities as cf) — s- 0. 

the Bethe-Yang equations have a solution corresponding to exceptional operators, while as 
far as <f> < (/> cr (A) the solution ceases to exist. A characteristic property of cf> cr (A) is that 
it vanishes in the limit A — > 0. Importantly, one finds that when <p approaches (f) CT from 
above the complex rapidities 1(2,3 move towards the branch points of the string tt-plane at 
— 2 =F i/g, where the function Y2 develops a double pole. On the z-torus the branch points 
correspond to the points of intersection of the boundaries of the string and (anti-)mirror 
region, see figure 1. Decreasing the twist below (fi CT , the only way to smoothly continue the 
evolution of «2 and U3 compatible with reality of the energy is to assume that they move 
along the cuts of the string u-plane or on the z-torus along the boundaries of the string 
region in opposite directions, reaching the positions of the exceptional rapidities at = 0. 
On the z-torus all the way towards the branch points the rapidities U2 and U3 remain 
complex conjugate but they loose this property upon passing them. On the u-plane this 
corresponds to the fact that 112 and 1*3 move along the lower edges of the cuts which reflects 
our choice of the string u-plane. In fact, such a behavior of 1^2,3 is the same as the one 
found in [11] for a two-particle BPS bound state at infinite J. Concerning the divergency 
of rapidities in the twisted theory, it is (almost) certain that this is just an artifact of the 
perturbative expansion. For finite A the rapidities may have an essential singularity at 
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(f> = <p cr such that the limit </>—>■ would produce the exceptional rapidities we conjecture. 
For example a term (f> e -A /(^ - ^ cr ( A )) 2 leads to poles in (f> in perturbative theory while for 
finite A it gives a zero contribution in both limits 4> — > 4> CI and <\> — > 0. It would be 
important to further justify the above-described scenario, in particular to construct the 
TBA equations for <\> < cr and show their consistency with our assumptions of positioning 
the rapidities on the boundaries of the string region. It is worth stressing that for these 
rapidities the usual asymptotic description does not exist because some S-matrices are 
singular. Nevertheless, the existence of the TBA for exceptional rapidities indicates that 
the corresponding construction must exist also for this case. 

The paper is organized as follows. In the next section we discuss the emergence of sin- 
gular states in the asymptotic Bethe Ansatz and introduce a twist. For the three-magnon 
case and L = 6 we also provide a perturbative solution of the Bethe- Yang equations up to 
the order A 6 accompanied by a small (^-expansion which reveals a singular nature of the 
state under consideration. In section 3 we discuss the twisted TBA equations for singular 
states and also compute the first Liischer correction to the energy for the state correspond- 
ing to Oq. We then show that the energy admits a smooth limit 4> — > 0. To shed further 
light of finiteness of energy in the twisted TBA approach, we explicitly demonstrate a can- 
cellation of the leading singularities in the expression for the energy at order A 7 . Section 
4 is devoted to the TBA approach based on exceptional rapidities. After formulating our 
conjecture on the exact form of Uj, we analyze the analytic properties of the asymptotic 
Y- functions which appear to be remarkably simple. Relying on the analytic structure of the 
asymptotic solution, we construct the corresponding TBA equations and show that they 
imply the fulfillment of the exact Bethe equations. We then compute the energy of Oq 
and show that in spite of the fact that in the approach based on the exceptional rapidities 
the TBA starts to contribute to the energy already at half- wrapping, the energy perfectly 
agrees with that found from the twisted TBA up to and including the first wrapping order. 
In the conclusions we discuss some interesting problems for future research. Some technical 
details are relegated to four appendices, and explicit expressions for twisted rapidities and 
Y-functions can be found in the Mathematica file attached to the arXiv submission of the 
paper. 

2 Bethe- Yang equations and singular rapidities 

In a perturbative expansion in g = t0 wrapping effects contribute to the scaling dimension 
starting from order g 2L where L is the length of the operator under consideration. Con- 
sequently, the Bethe- Yang equations provide the description of the perturbative spectrum 
up to the first wrapping order, and its predictions are usually expected to be qualitatively 
true even for finite but small g. It is therefore natural to start our analysis of exceptional 
operators with the corresponding Bethe- Yang equations. 

In what follows we will interchangeably use the gauge and string theory language, 
speaking equivalently of scaling dimension (of a gauge invariant operator) and energy (of 
the correspondent string excitation), etc. 
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2.1 Singular rapidities in the one-loop Bethe Ansatz 

The one-loop spectrum of M = 4 SYM in the su(2) sector is described by the XXX spin 
chain [3]. Scaling dimensions can be found by solving the Bethe ansatz equations for 
rapidities of M magnons 

Invariance under cyclic permutations 5 implies 

M . M 

e iP = l TT l±2^ = i W ith P = y p ( Uk ), p(u) = -i log — . (2.2) 

fe =i ? - Ufc fa 

The one-loop scaling dimensions, or energies, are then given by 

A/ 



fc=i 

Solutions of the Bethe- Yang equations exist also for complex values of the rapidities. It 
has been observed [4, 5] that among those there exist solutions with odd M where three 
rapidities are placed at 

ui = 0, U2 = — i , U3 = i, (2.4) 

and the remaining M — 3 rapidities come in pairs. The first three rapidities are rather 
exceptional: the corresponding momenta read 

7T 7T 
Pl=7T, p 2 = --+tOO, p 3 = ---lOO, (2.5) 

and similarly the individual energy of each of the last two magnons is ill-defined, signaling 
the necessity to introduce a regularization. This can equivalently be done by introducing 
a regularization parameter e in the solutions u\ = /i(e) and 112,3 = =R + /2,3(e) as in [2, 5] 
or by introducing a twist (j) in the Bethe- Yang equations as e.g. in [9]: 

l = e-«(^) L f[ Uk k = l,...,M. (2.6) 

Then the cyclicity condition (2.2) becomes P = M(j)/L mod 2ir. 

Focusing on the case M = 3, where only the three exceptional rapidities are present, 
one finds that when L is even (and of course L > 6) solutions can be found so that in the 
limit (f> — y rapidities tend to ui = and 112,3 = T*- This can be done by requiring that 
the divergence of momenta for small (ft is compensated by a singularity in the S-matrix 
Sxxx(u2, U3). Schematically one then has 

ui ~ (ft , U2 ~ — i — (ft — i 4> L , U3 ~ +i — (ft + i (ft L , (2-7) 



3 In string theory this is equivalent to imposing the level-matching condition. 
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where the value as well as the sign of the coefficient of the imaginary correction depends on 
L. Then, for all even L, the scaling dimension of the operator or equivalently the energy 
of the dual string state is also regular and reads 

lim E(4>) = L + 3g 2 . (2.8) 

Furthermore, the corresponding one- loop eigenvectors of the dilatation operator can be 
found by taking the limit of the Bethe wave-function of the twisted solution, yielding the 
N = 4 SYM operators (1.1). Therefore, at one loop, we conclude that there exists a 
family of eigenstates of the dilatation operator that can be constructed out of a building 
block of three exceptional magnons. These can be thought of as one magnon of maximal 
momentum p± = ir and one "infinitely tight" two-magnon bound state having maximal 
momentum P2 + P3 = — vr. It is interesting to see whether and how this picture changes 
beyond one- loop. 

2.2 All-loop Bethe- Yang equations and their breakdown 

The all- loop Bethe- Yang equations in the su(2) sector [2, 47, 48] including the twist 6 read 

M 

where L = J + M and cr(u k ,Uj) is the dressing factor. Here and in what follows we adopt 
the notation usual to field theory in which rapidities approach constant values for small g. 
Therefore, 

x± = x s {u k /g ± i/g) , x s (u) = |(l + sfl - 4/n 2 ) , (2.10) 

and the relation between rapidity and momentum of a magnon is e iPk = x^/x^. Again, 
the equations are supplemented by the level- matching condition e tP = z lM< &l L . 

As before, we focus on three-excitation solutions that for small g tend to the one- loop 
configuration of the previous section. From field theory, one expects the scaling dimension 
of any operator to admit a well-behaved small coupling expansion. Therefore, one would 
hope to resolve any singularity in the Bethe ansatz description by the same means used in 
the previous section. 

Let us consider, for simplicity, the case of the shortest operator of length L = 6. Then, 
for any non- vanishing value of </>, we can numerically solve (2.9). Some of these solutions 
are plotted in figure 2. These describe one particle with real rapidity and a pair of particles 
with complex conjugate rapidities for small g. However, as noticed in similar cases [27], 
it appears that the solution predicted by the Bethe- Yang equations breaks down at some 
critical value of the coupling g cr ((j)), which depends on the twist, see figure 2. There the 
rapidities are no longer complex-conjugate to each other, and as a result the energy becomes 
complex. 

6 As discussed in more detail in appendix 6.1, the twisted Bethe- Yang equations (together with the 
twisted level-matching condition) describe a 7-deformation of M = 4 SYM. 



(2.9) 
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Figure 2. Plots of the real and imaginary parts of M3 as functions of g for various values of <j>. For 
any fixed <p the rapidity U3 moves to the branch point — 2g + i in the field theory normalization 
and reaches it at g CI (cj>). Inset represents the inverse function 4> CT (g) which apparently is a linear 
function of g with slope ~ 3. 



We expect the breakdown to be an artifact of the asymptotic nature of the Bethe-Yang 
equations. What is striking, and peculiar of these states, is that the value of g C r{4>) where 
the breakdown happens goes to zero with eft, and therefore for finite g the twist cannot be 
removed no matter how small g is. This scenario also holds for larger values of L. 

This raises the question of whether the asymptotic description can be employed at 
least perturbatively in g. Expanding (2.9) perturbatively, up to the order g 2L ~ 2 one can 
find a solution of the form 7 



L-l 

11 

n=0 



l = ^/ l , n (0,L) 5 2 " + O( 5 2L ), (2.11) 



where at = the coefficients /j )Tl (0, L) are regular. The energy up to 0{g 2L ) is then 
found from the asymptotic formula 



M 



E ,s ym = j + Y j ^Jl + 4g2 sm 2 (p fc /2) , (2.12) 
k=l 

which involves the all-loop dispersion relation only. On general grounds we expect the 
asymptotic formula to receive corrections at order g 2L due to wrapping effects, and therefore 
to differ from the "true" result (which in principle might be computed by field theory 
perturbative techniques). For these particular states, however, the asymptotic energies 
appear to be divergent in the limit <f> — > at the wrapping order. For instance, in the 
L = 6 case we find 

£a Sym = 6 + 3,2 _ ^4 + 63^6 _ ^l g s _ _9_ (8C(3) _ 78% io + (213) 

/ 2187 3645 189783 81 27 1223982387\ 12 14 

+ V 1024^6 " 8192 04 + 131O72O0 2 + 128 C( ' + 32 C( ' ~ 14680064 J 9 + ^ ' ^ 



7 The solution for rapidities for L — 6 can be found in the Mathematica file attached to the arXiv 
submission of this paper. 
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Starting from the wrapping order g , the rapidities also become divergent in the limit 
0—t-O. This result is remarkable. Indeed, doing perturbative computations in 7-deformed 
N = 4 SYM one would find that for small (ft the numerical discrepancy between the 
asymptotic prediction and the true result is enormous. Obviously this is related to the 
fact that wrapping corrections have been neglected so far. Since the asymptotic energy 
diverges as (ft approaches zero, contribution of wrapping diagrams becomes crucial for 
diagonalization of the mixing matrix. This means that for exceptional states (or for states 
containing the three exceptional rapidities) a separation of the exact energy into asymptotic 
and wrapping parts is ill-defined in the limit of vanishing twist. 

In order to properly account for wrapping effects, we will use the mirror TBA. A 
convenient approach to excited states TBA is to make use of the contour deformation trick 
and of the knowledge of analytic properties of asymptotic Y-functions. For this purpose 
it is convenient to formulate TBA equations in the twisted theory for g < (ft where the 
asymptotic description can be trusted. 

3 Twisted TBA 

We want to find the mirror TBA description of the exceptional three-magnon configurations 
discussed in the previous section, which we expect to exist for any even L > 6. Our strategy 
will be to introduce a twist (ft and first formulate the TBA equations for the twisted theory, 
which corresponds to a 7-deformation of M = 4 SYM. 

Fixing a length L, for any nonzero (ft and for g small enough we can find the asymptotic 
solution of the twisted Bethe-Yang equations (2.9). These in turn allow one to write down 
the asymptotic Y-functions in the twisted theory. The details of this construction are 
given in appendices 6.2 and 6.3. Knowing the analytic properties of the asymptotic Y- 
functions, we can write down the TBA equations, which can then be solved numerically or 
perturbatively in g. 

3.1 Analytic structure of Y-functions 

We are considering here a family of configurations (labeled by even L) with one real ra- 
pidity ui and two complex-conjugate ^2,3, depending on g and (ft. Since eventually we 
are interested in the limit (ft — )■ 0, we restrict ourselves to considering a small region of 
parameter space, 

9 < eft « 1 , (3.1) 

where the first inequality follows from the necessity of having a real energy solution of the 
Bethe-Yang equations. 

Different states in the family have slightly different analytic structure for auxiliary 
Y-functions, that in turn yield different driving terms in the TBA equations by contour 
deformation trick. The procedure to formulate these equations in the case of complex 
rapidities has been detailed in [27], and can be applied straightforwardly to our case with 
minor L-dependent modifications. 
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Therefore, rather than attempting to give a unified description of each state in the 
family, we focus on the shortest one, with L = 6. In order not to clutter our treatment 
with technicalities, we relegate the discussion of roots of auxiliary Y-functions and the 
formulation of the TBA and exact Bethe equations to appendices 6.2 and 6.3. There we 
also briefly comment on how to obtain the TBA system for L > 8. 

Here, instead, we focus our attention on some peculiar properties of Yn functions for 
states with complex rapidities, which were also found in [27]. A crucial observation there 
is that depending on the location of the rapidities on the z-torus some Yg-functions may 
have poles inside the analyticity strip. As a result, there is a root of 1 + Yq located in 
the vicinity of a pole. If the rapidities lie just outside the analyticity strip, this leads to 
the appearance of extra terms in the TBA equations as well as the dispersion relation and 
total momentum quantization condition. 

This is precisely what happens in the case L = 6 for Yi . Let us indicate from now on the 
rapidities of the magnons as up . They obey the exact Bethe equations 1 + Y\*(up) = 0. 
Since we have for Y% that 

y 2 (4 1)+ ) = oo , Y 2 (uf') = oo , (3.2) 

and 4^ + and 4 are close to the real line then there exist two complex conjugate roots 
close to u^l such that 

l + y 2 (4 2)+ ) = 1 + Y 2 (4 2) ~) = 0. (3.3) 

Similar relations can be written also for Y3 close to 4^ ++ , but as it turns out, in the case 
of rapidities just outside the physical strip we can cast the TBA equations in a form that 
depends only on the usual roots and the (shifted) roots u^\- 

Taking e.g. the first equality in (3.3) and expanding Y 2 (u) = Res Y ?y£ around the pole 

u—u 2 

at 4 , one gets 



- (4 2) - u?) = Res Y 2 (4 1} ) + (4 2) " 4 1} ) + • • • • (3-4) 

For small residue of Y 2 this relation implies that — is of order of Res Y 2 which for 
small g is g 2L . It is also worth noticing that due to the presence of the poles (3.2) which 
are very close to the real line and almost pinch it, Y 2 (u) will take large values around 
u = Re(n2 1 ^) . 

3.2 Wrapping corrections for L = 6 at 0{g 12 ) 

We are interested in the first correction to the energy, which can be found from a pertur- 
bative expansion of the energy formula [27] 



j+t^)-hf: C dud t log{i+YQ) (3 - 5) 



3 

E 

i=l "" Q=l ' 

-^(4 1)+ )+^(4 2)+ )-^(4 2) -)+^(4 1) -) J 
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where we used the fact that for L = 6 the rapidities lie just outside the analyticity strip. 
To compute £(u^) to the order g 12 , it is sufficient to consider the asymptotic expression 
of the rapidities found by solving (2.9) and one obviously reproduces (2.13) from the first 
two terms in (3.5) since they correspond to (2.12). 

The leading perturbative correction due to wrapping effects can be found by expanding 
the remaining terms, 

= -±±£*%Y$ (3.6) 



Q=l 



-Z— — U. 



)Resy 2 °(4 iJ+ ) +z^p(n^~)Resy 2 °(4 ij_ ) , 



where we made use of (3.4) and replaced everywhere Yq by its asymptotic expression Yq, 
which can be found in appendix 6.2. Furthermore, at this order only the one-loop rapidities 
uf^ are needed. 

The final result is similar to the correction one would naively expect from Liischer's 
formula, with the important addition of the terms in the second line which are dictated by 
the contour deformation trick. It is worth noticing that, since Yq(u) > 0, the contribution 
of the first line alone is negative and for this reason can never cancel the small <f> divergence 
in (2.13). 

In the case L = 6 the computation of Ai?( wrap ) can be readily performed. As discussed 
above, the separation between the poles of y 2 ° a ^ ^4 ano - u 3 vanishes as <ft G for small 
g, as indicated by (2.7). Thus, the contributions divergent in the limit cj) — > come from 
the integral of Y" 2 ° an d from the residues on the second line of (3.6). Computing A.E( wrap ) 
and adding it to the asymptotic contribution, one finds that all divergent terms cancel out, 
giving in the limit <f> — >• the following result 

E = 6 + 3, 2 - \g* + - |V - JL(8C(3) - 78% 10 (3.7) 

/ 567 ., n , 189 .. r . 243 .,„, 84753\ 12 14 
+ (" 128 C(9) + 1M C(5) + 128 C(3) " 1024-J 9 + ° {d *> " 

(2) 

The cancellation of the divergencies would not be possible without the terms involving n 2 g. 
This provides the first justification of the energy formula (3.5) which does not rely on the 
contour deformation trick. 

3.3 Comments on the g 14 correction 

The cancellation of the divergencies at g 12 indicates that, when wrapping corrections are 
properly accounted for, the energy should not suffer from any singularity even at higher 
loop orders. On the other hand, considering the solution of the asymptotic Bethe ansatz 
(2.9), we find that not only the energy at g but also the rapidities at g 12 are divergent 
when the twist is removed. The mirror TBA is expected to render at least the energy 
formula finite. 

Unfortunately, even for the simplest L = 6 state, computing exactly the wrapping 
correction to the energy at order g 14 is a non-trivial task, conceptually similar to finding 
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the five- loop energy of the Konishi multiplet [30, 31], but much more involved because of 
the sophisticated analytic structure of the TBA system under consideration. 

To progress with the calculation of the energy at g 14 , one needs to know the rapidities 
uf^ at six loops. These cannot be found just by solving the Bethe-Yang equations: one 
has to consider the exact Bethe equations 

logY u (uf ) ) = (2n + l)Tri, neZ. (3.8) 

These equations are spelled out explicitly in appendix 6.3 and they involve auxiliary Y- 
functions as well as their roots. In a perturbative expansion, the exact Bethe equations 
can be written as 

= logFi*^) - (2n + l)7r* = logBY^u^) + 511^ (u^) , (3.9) 

where BY^ represents the Bethe-Yang contribution for particle i and 51Z^ is a correction 
of order g 12 (which also depends on the other rapidities, auxiliary Y- functions and roots). 
If we expand the rapidities u- 1 ^ around the asymptotic solution u°, 

u^=u° + 6uf-\ (3.10) 
we find that the exact Bethe equations can be rewritten as 

= E Su k ] + «* (0 («i) + 0(g u ) , (3.11) 

k=i k 

where we used that by construction BY®(u°) = 1. 

These three coupled equations are supplemented by the quantization condition of the 
total momentum P = 2irm + M<j>/L, where the total momentum is given by 

- i£ 2 (4 1)+ ) + i£ 2 (4 2)+ ) - ^ 2 (4 2) ~) + *£ 2 (4 1)_ ) • ( 3 - 12 ) 

Notice that the quantization condition is non-trivial because unlike most other cases, e.g. 
that of the Konishi operator where the rapidities come in pairs of opposite signs, P cannot 
be immediately seen to vanish due to the parity properties of Ig-functions. 

A natural question one may ask is whether the wrapping corrections to rapidities 6u\^ 
eliminate the divergent contributions in the asymptotic result at g 12 . In that case, it should 
be 

5uf ] = - (divergent part of u°) + 0(<jP) . (3.13) 

Without having to solve the complicated set of equations (3.11), we can plug our guess 
(3.13) into the total momentum quantization condition and check whether it is satisfied. 
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The advantage of this strategy is that at the order g 12 we can expand (3.12) as 



- i^( u °2 + ) Res Y 2 (u°+) + Res Y 2 (u°-) + 0(g u ) . (3.14) 

where the only non-asymptotic objects appearing in (3.12) are precisely o~ui\ 

Surprisingly, we find that the guess (3.13) is incompatible with the total momentum 
condition; in fact it would make P divergent as (ft 0. This implies that the individual 
rapidities found from exact Bethe equations remain divergent in perturbative theory. 

The only way of checking whether the g 14 wrapping correction to the energy makes it 
finite for small (ft is to deal with the full set of TBA equations and expand them around 
the asymptotic solution and then in powers of g. This is straightforward but cumbersome, 
and is done in appendix 6.4 for the case L = 6. The linearized TBA system ends up to be 
more complicated than in the case of the Konishi operator. In particular, the linearized 
system for the correction to Yj^t vw -functions does not decouple from the other auxiliary 
equations, which makes it hard to find an analytic solution. 

On the other hand, if we focus on the most ^-divergent part of the corrections to 
rapidities (which in turn determine the most divergent part of the corrections to the energy) 
it is relatively easy to see that once again the wrapping effects precisely cancel the asymptotic 
divergence. The compatibility of this cancellation with (3.12) can also be seen as a non- 
trivial check of the formula for the total momentum. 

In conclusion, we find strong evidence of a general mechanism by which the TBA 
description of the exceptional operators can be obtained by introducing a twist eft as a 
regulator. Even if the TBA system can be found from the asymptotic data only when 
g < (ft, and therefore never, strictly speaking, at (ft = 0, the resulting physical predictions 
will be regular in (ft when wrapping effects are accounted for. Therefore, we can compute 
the perturbative energy for small (ft and then take the limit (ft —> in the final result. 

Even if in principle a similar strategy could be repeated to find energies at finite g, this 
would require to (numerically) solve the full TBA system for several values of (ft in order 
to extrapolate to (ft — ¥ result. This would be practically unfeasible, and it is therefore 
important to look for an alternative TBA description of these operators, which does not 
resort to introducing a regulator. 

4 TBA with exceptional rapidities 

The twisted TBA approach provides a way to compute the anomalous dimensions of ex- 
ceptional operators in perturbative gauge theory. However, it leaves open a question of 
determining the dimensions at any value of the coupling constant. In this section we pro- 
pose a set of TBA equations which allows one to calculate the dimensions of these operators 
at any value of A. 

The main idea is that since an exceptional operator is dual to a string theory state 
which is composed of a fundamental particle and a two-particle bound state with maximum 
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i(Q-2)/g, -i(Q-2)/g 



Table 1. Relevant roots and poles of asymptotic Y- functions within the mirror region. 2 means 
either a double zero or a double pole at 0. 

allowed momenta ±7r, the Bethe roots in the gauge theory normalization for any exceptional 
state are in fact independent of the coupling constant: u\ = , u 2 = — i , u% = i. The roots 
«2,3 satisfy the bound state condition, and since their real part is 0, they are on the cuts of 
functions. According to [11], they must lie on the same sides of the cuts, and therefore, 
we propose that the exact Bethe rapidities (in the string theory normalization which will 
be convenient to write the TBA equations in this section) for any exceptional state are 
equal to 

ui = +— , u 2 = iO, u 3 = iO. (4.1) 

2 5 9 

With this choice of the signs in front of iO, the fundamental particle and the bound state 
composed of ^2,3 have momenta +7r and —tt respectively, if one uses Mathematica's con- 
ventions for branch cuts. Then the root u 2 lies in the intersection of the mirror and string 
regions, and U3 is in the intersection of the string and the second mirror regions. Notice 
that it is different from the state analyzed in [27] where the rapidity M3 was in the intersec- 
tion of the string and the anti-mirror regions. The location of the rapidities on the z-torus 
is shown on figure 1 , and in terms of the z-rapidity variable all Y- functions and dispersion 
relations are meromorphic in the vicinities of these points. 

These rapidities lead to a quite simple analytic structure of asymptotic Y-functions 
with double poles and zeroes at the origin of the mirror n-plane, see Table 1, and it is 
natural to assume that the exact Y-functions would have the same analytic properties. 8 

4.1 TBA equations 

In this subsection we list the simplified and hybrid TBA equations for the exceptional 
states. They can be obtained from the ones discussed in [27] by sending the roots tm 

8 Let us mention that Y-functions with double poles and zeroes at the origin are typical for boundary 
TBA, see e.g. [49-51]. 
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to 0, and uf' ,uf^ to itj. The only exception is the hybrid equations for Yq where one 
should take care of the fact that the root U3 is located in the intersection of the string 
region and the second mirror region but not in the anti-mirror region as it was in [27]. The 
TBA equations below are consistent with the analytic structure of Y-functions in Table 1 
supplemented by the conditions Yi„(0) = Yi t (—i/g) = Y\^(i/g) = —1. 



Simplified equations for Ym\ 

■ * , \ , 1 ft 1 T 7" \ f t 1 ~\ 7~ \ C , L ~ Y. 

j 

y+ 

Simplified equations for Y M \ VW 



logYjtflto = 21ogS'(-+v)+log(l + y M _ 1 | u ,)(l + y M+1 | u ,)*s + (5A/i log - - -y- *s. (4.2) 



loglAflDto =2hn logS(- + v) + log(l + Y M -l\vw)0- + ^Af+llw) * s 

1 - Y- 

+ S M1 ]og- — *s-log(l + y M +i)*s- (4.3) 

1 — X + 

Simplified equations for Y± 

log -± = log(l + Yq) -k K Qy - log S Uy ( Ui , v) , (4.4) 



Y. 

lo. 



g y + y_ = 2 log 1 * s - log (1 + y q ) *k q + 2 io g (i + y q ) * iffi * 8 

1 + Y l\w 

- 41og 5(i + v) - log Sl f 1(u ^ v)2 * s . (4.5) 
g ^ S 2 {ui-v) 

It is worth mentioning that since the driving terms in the equations above satisfy the 
discrete Laplace equation 

1 1 

S Q (v - -)S Q (v + -) = Sq-i(v)S q+1 (v) , S (v) = 1 , 
they can be written as 

- J2 lQ g Su v (ui,v) = - log S Uy (0, v) - log S 2 . v (0, v) , (4.6) 

i 

- E lQ g %7 { * * = ~ log - J , * a - 2 log (0, P ) * S + log S 2 (0 - v) . 

I 2 \Ui-v) 52(0 -V) 

This shows that the driving terms in eqs. (4.4,4.5) can be understood as appearing not 
because of the zeroes of 1 + Yi, at u = , — i/g , i/g in the string u-plane but due to the 
zeroes of 1 + Y\ t and 1 + Y 2t at u = in the string u-plane. It is consistent with the 
interpretation of an exceptional state as a bound state of a fundamental particle and a 
two-particle bound state with rapidities equal to 0. This interpretation however requires 
using integration contours different from the ones described in [27]. 
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Simplified TBA equations for Yq 



• Q > 3 

log yq = n S vT7^ * s • (47) 

• Q = 2 

(! + fH 2 

log F 2 = -2 log S( l n ~ v) + log n \' l ; w /^ *p.v s ■ (4.8) 
9 I 1 + y 1 Jv 1 + y 3 J 

Hybrid TBA equations for Yq 

To make the presentation transparent, we introduce a function which combines the terms 
on the right hand side of the hybrid ground state TBA equation (Ltba = J + 2) 

G Q {v) = -L TBA £ Q + log(l + Y Q ,)*(K^ ) +2s*KZx 1 ' Q ) (4.9) 
+ 2\og(l + Y llvw )*s*K y Q + 2log(l + Y Q _ llvw )*s 

~ 2 lQ g * s * + lQ g rTTT + lQ g (! " y) (! " ^) *^s/Q • 

With the help of Gq , the hybrid TBA equations for Yq read as 

\ogY Q {v) = G Q (v) -Y,logS 1 s ? { Q ) (u i ,v)+4\ogS* p . v . K% x (--,v) (4.10) 

i i 2i 2i 

- logS Q (-^ - v)S y Q(-^,v)S Q (-v)S yQ (0, v)Sq( — - v)S yQ ( — ,v) . 

It is important to stress that since the location of the Bethe rapidities is exactly known 
the only parameters in the TBA equations for exceptional operators are the charge J (or 
equivalently the operator length L = J + 3) and the coupling constant g. In this respect 
these TBA equations are of the same level of complexity as the ones for the ground state 
of any integrable model. 

4.2 Exact Bethe equations 

To construct the TBA equations by using the contour deformation trick one has to assume 
that 1 + Y± t has zeroes at u = ,—i/g ,i/g in the string plane. On the other hand once 
the equations have been derived one can use the analytic continuation to calculate Y± at 
these points. Thus, the conditions 

Fi.(0) = -1, Y u (--) = -l, Y u (-) = -l, (4.11) 
9 9 

on Y\ t must follow from the TBA equations. This imposes nontrivial consistency conditions 
on the TBA equations which we discuss in this subsection. 
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Bethe equation at u% = 0: 11,(0) = —1 

We begin by showing that 11,(0) = — 1. Indeed analytically continuing the equation for li 
to real v one gets 

log Y U (v) = G u (v) - ^V ( 2) K v ) 

i 

+ 41og Res 5 * Klt x (--,v) + 2 log S^(--,v) - ilog(-v - - ^l^ % W 

9 9 9 

i i 2% 2% 

- log5i(-- - v)S yU (--,v)S 1 (-v)Syi,(0,v)S 1 (— - v)S yU ( — ,v) . 

Then one finds that the imaginary part of Gi t (v) in the limit v — > is equal to itt(J + 2) 
because all the kernels in Gi t (v) are antisymmetric at v = 0, and the real part of G\ f (v) 
is given by the usual expression 



Re G u (v) = -J2 log S Uv (ui, v)*Ki. (4.12) 



One can then easily check that in the limit v — > 

i i 

i i 2i Xs ~( \ 

+ 41ogRes S*K% X (— ,v) + 21og^(— ,v) - 41og(-t, - -) _ % W 

9 9 9 x °W~^k) 

i i 2i 2i 

-log Si ( v)S v u( — ,v)S 1 (-v)S yU (0,v)S 1 ( v)S y i t ( — ,v) = mod 2m , 

9 9 9 9 

and therefore 

logY u (0) = ivr(J + 2). (4.13) 
Thus if J is odd as it is for exceptional operators then ii*(0) = — 1. 
Bethe equation at U2 = —i/g'- Y\{—i/g) = —1 

To show that Yi f (—i/g) = —1 we notice that U2 = —i/g — i0 is in the mirror region, and 
therefore li„ (112) = Y\(v,2)- Moreover, since we approach —i/g from the mirror real line, 
we can always use the mirror-mirror kernels in (4.10). Then to show that Y\(—i/g) = — 1 
we use that all Y- functions are even, and all the kernels in (4.10) satisfy 

K(t, v) = K(-t, -v) , (4.14) 

and therefore for any even function / 

2/ * K(v) = f * K(v) + f * K(-v) = / * (K(v) + K(-v)) . (4.15) 
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Thus we have the following equality 

21o g yi(u) = G x {v) + Gt(-v) - 2^1og5 s 1 ( ^ ) (n i ,«) (4.16) 



+ 41ogS* p .„. (Kll x (--,v) + K%, (-1 - v )) 

9 9 

i i 2i 2i 

- 21og5i( v)S y i( — , v)S 1 (-v)S yl (0,v)S 1 ( v)S yl ( — ,v) . 

9 9 9 9 

Now we want to take the limit v — > —i/g. Since all the kernels satisfy the discrete Laplace 
equation we would naively get 

di-i/g) + G^i/g) = G 2 (0) - log (1 + Y 2 (0)) , (4.17) 

where the last term appears because of the pole in K 2 ^ 2 Jt,v) at t = ±i/g. The kernel K y Q 
also has a pole there and it produces the term 2 log (l + Yi\ vw ) * s which is in G2, and it 
could produce the term log (l — y— ) (l — y— ) but it vanishes because Y±(0) = —00. The 
only problem with (4.17) is that 12(0) = 00, and therefore we should deal with the term 
F\ = — log (1 + Y2) * K\ more carefully. We represent it in the form 

F x {v) = - J dt\og 1 ^ 2 ^ K^t-v)- f dt log (l + ^)Kx(t-v), (4.18) 
where C 2 = lim^o t 2 Y 2 (t). The first term then represents no problem and one gets 



2F 1 {e --)=-/ dt log (1 + Y 2 {t))K 2 {t - e) (4.19) 



/C 2 f C 2 

dt log (1 + -p)K 2 {t - e) - 2 / dt log (1 + -p)K x (t - e + 



9 



where e is infinitesimally close to with positive imaginary part. The integral on the 
second line can be computed, and expanding it in powers of e one gets 

i C 2 
2Fi(e - -) = - log (1 + Y 2 ) * K 2 (0) - log - in . (4.20) 
9 e 

Thus, the formula (4.17) contains the extra in term, and takes the form 

Gxi-i/g) + Gx(i/g) = G 2 (e) - log Y 2 (e) - in + o(e) . (4.21) 
Taking into account the TBA equation for Y 2 one gets 

21ogYi(w) = -itt + £ log («<,« + -) - 2 ^ log S^ (^) (4.22) 

i 9 i 

i i 2i 2i 
- 2 log Si ( v)S y i( — ,v)Si(—v)S y i(0,v)S 1 ( v)S y i(—,v) 

i i 2% 2i 

+ logS 2 (-- - v)Sy 2 (--,v)S 2 (-v)Sy 2 (0,v)S 2 (— - v)S y2 ( — ,v) , 
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where v = e — ^. Taking the limit e — > one finally gets 

logYi(--) = -in. (4.23) 
9 

In the same way one can show that Yi(^) = — 1 (or one can use the Y-system equation 
for Y±), and then the condition Yi»(g) = — 1 can be proven by using the crossing symmetry 
relations as was done in [27]. Let us finally mention that it should be possible to show that 
the TBA equations imply in addition ^2,(0) = — 1 because the particles with rapidities 
±i/g can be thought of as constituents of a two-particle bound state with rapidity equal 
to 0. This however requires a careful analytic continuation of the hybrid TBA equation for 
Y2 to the string u-plane through the cut at —2i/g, and we will not pursue this here. 

Scaling dimensions of exceptional operators 

Scaling dimensions of exceptional operators or energies of dual string states are found from 
the usual formula 

A-J = E- J= J°° dn^log(l + y Q ) 

= yi+V + + V - ^ I" du ^ Ml + Yq) , (4.24) 

where we used the exceptional rapidities of the particles. This formula shows that at large 
g the first two terms in (4.24) which come from the dispersion relation are proportional 
to g. On the other hand for finite J and large g the scaling dimension of these operators 
should behave as yfg. Thus, the linear term should be canceled by the contribution coming 
from the Ig-functions. This is different from the expected large g behaviour of two-particle 
states studied in [25, 29]. It would be interesting to understand if the linear term comes 
entirely from the pole contribution of Y%. 

4.3 Leading TBA correction up to g 10 

The proposed TBA equations are based on the assumption that the rapidities of exceptional 
states are given exactly by (4.1). These rapidities are obviously very different from the 
rapidities of the states in the twisted theory which diverge in the limit </>—>■ at least in 
the perturbation theory. Still, the TBA equations should produce the same perturbative 
expansion of the scaling dimensions of exceptional operators as the one we obtained from 
the twisted TBA equations in the previous section. In this and next subsections we compute 
the scaling dimension of the shortest exceptional operator of length L = 6 and show that 
it coincides with the twisted TBA result. We will use the gauge theory normalization of 
rapidities in which the exact Bethe roots are 0, ±i. 

Let us recall that the finite-size corrections to the energy of the twisted exceptional 
operator for finite (f) start exactly at g 12 as expected for an operator of length L = 6 from 
the 5u(2) sector. Thus up to g 10 one can just use the dispersion relation and the BY 
equations. Then, as was shown in the previous section, one gets 

j? 9 5 4 63g e 621 5 8 9g 10 C(3) 7047g 10 

^= = 6 + 35 + 32— + ^56~- (425) 
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On the other hand if one uses the energy formula (4.24) with the exceptional Bethe roots, 
then the contribution coming from the dispersion relation is just given by the first two 
terms and its expansion up to g 10 produces 

pasvm n z-R r. y-j n „ 2 9g 4 33g 6 645g 8 3591 5 10 tA ne . 

£.asym = ^/\ + 4„2 + ^4 + 4gf 2 ^ g + 3^ ^_ _| » ^_ _| ^_ (4.26) 

v y v y y 4 8 64 128 v 1 

The two formulas obviously become different already at the order. Thus the finite-size 
corrections in the case of the TBA with exceptional rapidities must appear at the g e order 
which from the field theory point of view is half-wrapping. We know that perturbative 
expansion of all Yg-functions begins at g 12 and therefore any YQ-function regular on the 
real line begins to contribute to the energy at the g 12 order. The only exception is Y 2 - 
function which has a double pole at zero (if eft = 0). As a result the perturbative expansion 
of the integral J du^ log(l + Y 2 ) starts at the g & order. Thus, up to the g 10 order one 
should get the same energy (4.25) by keeping only Y% in TBA equations and the energy 
formula. Therefore, the formula of interest up to g 10 is 

E = £asym _ _L f d J^ j (1 + y 2 ) (4.27) 
Z7T J au 

where £J asym is given by (4.26). Up to the g 10 order we only need the coefficient of the 
double pole at u = up to the g ie order 

v( , 9g 12 (3g 4 (8C(3) + 15) - 24g 2 + 8) , + n( 2, 

Y 2 (u) = + const + 0{u). 

Then computing the integral in (4.27) one finds 

E P* = _ J_ / „ „* log(1 + y 2) = _ V°«3) _ 135^ + V _ 31 

2vr J dv BV ' 32 256 8 16 v ' 

where in Y% we only kept the 1/u 2 term. 

Adding (4.28) to (4.26), one gets precisely (4.25). 

4.4 Next-to-leading TBA correction at g 12 

The agreement between the energies observed in the previous subsection should also hold 
at the g 12 order where one should calculate the usual contributions from all Yg-functions. 
In addition one also has to take into account the TBA correction to the coefficient of the 
double pole of Y 2 which is of the g 18 order. 

Linearization of the TBA equations 

It is well-known that at small g Y-functions get TBA corrections beyond their asymp- 
totic form Y°. Computing the leading TBA corrections requires linearization of the TBA 
equations which can be done by representing any Y-function as follows 

Y(u) = Y°{u) (l + &{u) \ . (4.29) 
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Since the Bethe roots do not get corrections, the $^'s have neither zeroes nor poles on 
the real line. Then one expands the hybrid TBA equations up to the first order in ^aux 
while keeping only the contributions from the asymptotic Ig-functions on the r.h.s. of the 
equations. It is clear that leading corrections to any & are of order g 6 or higher, and they 
come only from the pole part of Y 2 °. Discarding any term of 0(g 8 ), we find that only the 
following two equations are relevant at the g 6 order 

^ 2 = log(l + Y 2 °) * (K% 2) + 2s* Kl 2 wx ) + 4 (A llvw & 1]vw ) * s , (4.30) 

★ s + A M+1 \ vw ty M+1 \ vw *s- S M i log(l + Y 2 °) * s , (4.31) 

where we defined the coefficient 

yo 

, M\vw j, , ^ 

A M\vw = -, V o ' M - l - 

1 + 1 M\vw 

The g 6 contribution of Y 2 ° to these equations can be easily computed because for any kernel 
K(u, v) regular for real u and v one gets 

log(l + Y 2 °) *K^R° 2 I du log(l + ^ ) K(R° 2 u, v) -> ^-g 6 K(0, v) , (4.32) 

where R 2 is the square root of the coefficient of the pole of Y 2 ° 

Y 2° = ^jjT + • • • . ^2 = i? 6 (l " + ^(8C(3) + 9) + ^(-120C(3) - 108C(5) - 55)) . 

(4.33) 

This also proves that the leading TBA corrections to Y- functions are of order g 6 . 
Expansion of the energy formula 

Let us now assume that we know $2 up to the g 6 order and compute the energy up to the 
g 12 order. The expansion of E asym gives 

£iasym, (12) _ 

512 

The contribution of Yq with Q / 2 is found from the usual formula 

Q^2 J 

Computing the integrals and taking the sum, one obtains 

E {Q+i) = 12 <^ 135C(3) 297C(5) _ 567C(9) 358424597369 \ 
y y V 128 128 128 580608000000 7 1 ' J 

To find the contribution of Y<i we represent the integrands in the energy formula as follows: 



1 + ^2 , - R 



log(l + Y 2 ) = log — ^ + log (1 + -f ) , (4.37) 



1 + 



„2 
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where i?2 is the square root of the coefficient of the pole of Y2 which also includes the 
contribution from an d therefore can be written as 

R 2 = RW 1 + ^2(0)- (4.38) 

The first term in (4.37) is regular everywhere, and can be expanded in g starting from g 12 , 
and at that order depends solely on asymptotic quantities. Its contribution to the energy 
at the g 12 order is given by 

rnreu 1 f , , 1 + ^2° 15795402631 12 . . 

Eg* = / du log -f- = g 12 . 4.39 

Y2 2ir J 1 + ^ 580608000000 J V ' 

The contribution of the second term yields 

< 6 = ~l / du f log (1 + §) = Eff - |!%(0) + Otf% (4.40) 
where E^t e is the contribution due to the pole of Y£ 

E$> = 3 ^(72 C (3) + 54 C (5) + 55) - 9 ^ - ™£ + ^ - 1. (4.41) 
y 2 256 v sv ; sv 1 ' 32 256 8 16 v ' 

This means that to find the energy at order g 12 , we need to know the leading TBA 
correction to Y% at u = 0. The correction is given by (4.30) which at u = can be written 
in the form 

3o 6 

^2(0) = ■^-(81og(2)-3) + 4(A 1 | tno 91| wo )* a (0). (4.42) 

The last term 4 (Ai\ vw * s(0) can be found by solving eq.(4.31) which takes the 

following explicit form 

3<? 6 

<&M\vw(u) = ^M-\\vv^M-\\ww * s + A M+l\vw%t+l\vw * s ~ Smi~^~^ s(u) . (4.43) 
Introducing the functions Xm{u) which satisfy the following difference equations 

— = %-i + %i + S M i 2-ks[u) , (4.44) 

the quantity 4 (Aiu w ( 3f\\ v J) * s(0) appearing in (4.42) can be written in the form 

3a 6 

4 (Aj^ S^) * s(0) = - X x (0) . (4.45) 

Thus summing up all the contributions one finds the energy of the exceptional state at the 
g 12 order 

(12) = 3g 12 (24JT 1 (0) - 1512C(9) + 1008C(5) + 648((3) - 28237 - 241og2) 

1024 ' 1 ' ' 
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Comparing this formula with (3.7) obtained from the twisted TBA, one gets 

E {12) - 4=o = ^<? 12 (12iTi(0) + 7 - 12 log 2) . (4.47) 
Thus the two results coincide if 

SC\ (0) = log 2 - — « 0. 109814 . (4.48) 

We could not prove this equality analytically. Solving the system (4.44) numerically we 
find that the equality (4.48) holds with very high precision. 

To conclude this section let us point out that the consideration above can be easily 
generalized to the exceptional operator of length L = J + 3. The Y^-function begins to 
contribute at the g L order. The improved dressing factor contribution can be easily found 
at this order, and one gets that the energy of the exceptional operator is just equal to 

E L = J + y/1 + 4 5 2 + y/l + V _ JL/ + 0{g L+2 ) . (4.49) 

It is not difficult to check that at this order the same expression is obtained by using the 
twisted state in the limit <p — > [2]. One can in principle go all the way till g 2L . The only 
technically nontrivial part is finding the power series expansion of the dressing phase up 
to the g L+2 order. 



5 Conclusions 

In this work we have provided the mirror TBA description for the exceptional class of gauge 
theory operators Ol- From the point of view of the Bethe Ansatz the states corresponding 
to these operators are singular that is the asymptotic energy diverges at the first wrapping 
order in the limit of vanishing twist. On the other hand, in the approach based on Baxter's 
Q-operator, the same state with M = 3 Bethe roots can be described by means of L — 
M + 1 = L — 2 dual roots which are all regular at one loop. It would be interesting to see 
whether the dual root picture can be implemented at the level of the TBA equations. A 
natural starting point here would be to explicitly develop the all-loop Baxter equation in 
the su(2) sector in the spirit of [52]. 

In a certain respect the operators from the family {Ol} are even more interesting 
than the Konishi operator. Indeed, the fact that their Bethe rapidities are known exactly 
must simplify the numerical analysis of the corresponding TBA equations since one does 
not need to solve the exact Bethe equations. Also, the rather rigid analytic structure of 
Y-functions - the presence of double poles and zeroes - hints that it possibly remains the 
same all the way from weak to strong coupling which might help to find a proper ansatz for 
Y-functions at strong coupling. This should be contrasted to the case of regular operators, 
where the position of zeroes and poles depends on the coupling constant and there are 
critical points [25, 29]. 
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Since a three-magnon state with rapidities 0, +i/g, —i/g can be viewed as a scattering 
state of a fundamental particle and a two-particle bound state with momenta ±7r, the 
asymptotic energy is 

^asym = J+ ^ + Ag2 + ^ 2 2 + 4 5 2 . 

Therefore, at large 3 the asymptotic energy scales as E RSyra ~ 5. On the other hand, the 
operators we consider belong to the class of short operators for which the energy must scale 
as y/9 ~ v^A at strong coupling. Hence, according to the TBA description, the contribution 
of Ig-functions must scale as g at strong coupling and cancel the leading term of i? asym 
at g — > 00. It would be interesting to verify this fact by constructing the corresponding 
analytic and numerical solution. 

Let us also mention that recently there has been an interesting development [53]- [55] 
concerning a construction of a finite set of non-linear integral equations (NLIE), which 
is a complementary approach to the TBA description of the spectrum of the AdSs x S 5 
superstring. It would be important to see how the states corresponding to operators Ol 
can be accommodated within the NLIE approach. 

The experience we gained here with the exceptional operators brings us back to the 
question of the strong coupling behavior of a generic bound state in M = 4 theory discussed 
in [27]. We expect that similarly to what happens in the <j> — > limit for twisted states, 
the complex rapidities of a generic bound state will reach the branch points at finite value 
of g and afterwards continue to move along the boundary of the string region towards the 
position of the exceptional rapidities reaching them at g = 00. To confirm this picture 
one has to further investigate the TBA equations obtained in [27]. If true this would 
suggest a universal behavior of a generic state: when coupling increases eventually real 
rapidities move towards —2,0,2, while complex rapidities reach the branch points and 
upon passing them approach the exceptional rapidities. The points —2, 2 and 0, ±~ would 
serve as attractors for all rapidities. This would classify states with a finite number of 
roots at strong coupling and might explain the universal v^A-behavior of the energy of 
short operators. 

Acknowledgements 

We are grateful to Niklas Beisert, Nadav Drukker, Gregory Korchemsky, Matthias Stau- 
dacher and Stijn van Tongeren for useful discussion. We also thank Stijn van Tongeren for 
useful comments on the manuscript. G.A. and A.S. acknowledge support by the Nether- 
lands Organization for Scientific Research (NWO) under the VICI grant 680-47-602. The 
work by G.A. is also a part of the ERC Advanced grant research programme No. 246974, 
"Super symmetry: a window to non-perturbative physics". The work of S.F. was supported 
in part by the Science Foundation Ireland under Grant 09/RFP/PHY2142 and by the 
Institute for Advanced Studies, Jerusalem, within the Research Group Integrability and 
Gauge/String Theory. 



- 25 - 



6 Appendices 

6.1 Twisted transfer matrices and relating twist to a 7-deformation 

In this section we will show how the twist parameter <p that we have introduced as a mere 
regulator can be related to the parameters of a 7-deformation of N = 4 SYM. To do this 
let us recall that the most general 7-deformation imposes twisted boundary conditions on 
the angles ipi of S" 5 as follows [8] 

99j(2vr) = ^j(O) - 2vr e ijk -/jJ k , (6.1) 

where 7^ are three deformation parameters, and J, are angular momenta on 5 5 correspond- 
ing to the direction of ipi . Let us introduce the notation 

a.i = -27T e ijk -/jJ k , a e = , a r = . (6.2) 

The level-matching condition in the presence of such modified boundary conditions is 

P = Qi + 2vrn, n£Z, (6.3) 

and the asymptotic su(2) transfer matrix in the left and right sectors have the form [43] 



M - - 1 / 7 + M _ + - 




i • I } ' ■ , ' '- 1 ■' . ■ j ' i ' i ' ~\~ ry* ry* I 

8=1 ii x + V * «=1 J i x+ 

where M is the number of magnons and x^ , xf are the usual parameterizations of mirror 
and string rapidities. 

We will restrict to the choice 

q 3 = 0, a = ae = a r = -—, (6.5) 

and it is immediate to obtain the Bethe-Yang equation 

-1 = W, (6.6) 

from the analytic continuation of the asymptotic Yq functions 

M 

Yq{v) = e- jiQ(v) (v, {u k }) T%1 („, {uk }) [J S^(v, Uj ) . (6.7) 

3=1 

One then finds 



M 



X k X j 



(6.8) 
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which can be rewritten using the explicit form of the S-matrix and the total momentum 
quantization condition (6.3) as 



e i Pk (J+M) e ia 2 e -<ai 77 u k ~ Uj ~ 2z 



M 



u k - Uj + 2i 

Applying this discussion to the family of the states of interest, for which M = 3, 
J = Ji = L — 3, J2 = 3 and J3 = 0, one finds 

3 7 i + (L- 3)72 = 0, (6.10) 

whereas the Bethe-Yang equations can be written simply as 

so that we can think of twist as being related to a deformation by 

= -2^L 73 = ^-^a. (6.12) 

It is also interesting to notice that, in the case L = 6, the constraint (6.10) is compatible 
with the choice 

71 = 72 = 73, (6.13) 

which is the Leigh-Strassler deformation preserving J\f = 1 supersymmetry and dual to 
the Lunin-Maldacena background [7]. Furthermore, inspecting (6.4) one finds that, on a 
solution of (6.9), the explicit dependence on the deformation parameter drops from the 
asymptotic transfer matrix. As a result, many of the analytic properties of the asymptotic 
Y-functions will be essentially the same as in the untwisted case. 

6.2 Twisted Y-functions and their analytic properties 

The asymptotic transfer matrices in the antisymmetric representation (6.4), together with 
Bazhanov-Reshetikhin formula [56], yield all of the Tqqi. 9 From those, one finds the 
auxiliary Y-functions [13] 

v Ti ) mTi > m+2 v r 2) i 22 )3 T 2) i Tm,iTm+2,1 f(i -lA\ 

*M\w — ™ > Y - — —7^ — > Y + — — 7f — i Y M\vw — — 7f j 10.14J 

J-2,M+1 J 1,2 -tl,2-i3,2 J-M+1,2 

whereas the asymptotic Yq functions are given by (6.7). All are real analytic functions of 
the mirror rapidity. The relevant analytic properties of the full Y-functions can be found 
from inspecting their asymptotic counterparts at small g. Recall that in doing so, we will 
always consider the regime cp < g. 

In table 6.2 the meromorphic structure of Y-functions is schematized. A few remarks 
on how this scenario depends on L are in order: 



9 For practical purposes it can be convenient to directly find T x qi by a duality transformation as detailed 
in [21] rather than from Bazhanov-Reshetikhin formula. 
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7pynpc 


Poles 


1 M\w 


r A M±l 

'Mil 




1 T" *M|lU 


A M - „A M + 
M ' ' M 


U2 — (ivi -\- 1)1/ g , 113 -f- (iw -f- L)i/ g 


y: . 
- 1 l\vw 


n l i r o 




1 _l_ V, ,, 
1 T i M\vw 




i /" /\/f i i\ n ' /a f ?\/f _i_ i \; //> 


V 
I 


+ 


U~2 , % 


y. 




"l 


1 - Y_ 


„Ao - „Ao + 
'0 ' '0 
















y 2 




n 2~ ' u 3 


*Q,Q>3 




u 2 +t(Q-l), U3 -i(Q-l) 



Table 2. Relevant roots and poles of asymptotic Y-functions for general L. The index Am = 
1, A™/" labels different roots, and A™ f ax depends on L. 

1. Auxiliary functions Ym\w an d Y_ satisfy quantization conditions at the shifted values 
of the (real) roots {r^j* Iam^,...^™") which by contour deformation trick will appear 
in the TBA equations. Their number A^ ax and their position will depend on the 
value of L under consideration. 

2. As discussed, the form of the TBA equation and of the energy and momentum for- 
mulae will depend on whether the complex rapidities 1*2,3 lie inside or outside the 
physical strip, which depends on L. 

3. As seen in the previous appendix, the case L = 6 is special in that it can be linked 
to a deformation which preserves more supersymmetry. As a result, the large- -u 
asymptotic of Yq(u) will be different depending on whether L = 6 or not, which is 
consistent with the fact that the relation between the TBA length Ltba and J is 
modified when all supersymmetry is broken [43]. 

4. It is worth pointing out that Y2 has poles at , which lie very close to the real 
line. As can be seen from (2.7), in the limit g <C <f> <C 1 their distance from the real 
line is of order <fi L . 

6.3 TBA equations for the twisted theory 

The TBA equations for the family of states of interest can be engineered by contour defor- 
mation trick, taking into account the analytic properties for the state at hand. We write 
them in a rather general form, by introducing terms D* that indicate the driving terms of 
a given equation that depend on the roots {r\f }\ M =i,...,\™ ax , coming from Y M \ W = — 1 or 
Y_ = 1. 

For concreteness, we consider a more involved case in which the complex rapidities 
lie (just) outside the analyticity strip (which is the case of L = 6), and express TBA 
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equation in terms of simplified and hybrid equations only. When the rapidities are inside 

the analyticity strip there is no need to consider the quantization of the roots of 1 + I2 

(2) 

and therefore u 2 3 drop out from all equations. We refer the reader to [27] for a detailed 
discussion of the TBA equations with complex rapidities, whereas the definition of the 
kernels used below can be found in [25]. 



Simplified equations for Yj 



M\w 



1-^ 



log Ym> = log(l + Y M _!\ W )(1 + Y M+1 \ W ) *s + S M x log 1 _ x * s + D M\-w ( 6 - 15 ) 
Simplified equations for Y M \ VW 



^ogY M \ vw = - log(l + Y M +i) * s + log(l + Y M _ llvw )(l + Y M+1 \ VW ) * s (6.16) 

S{u { 2 2)+ -v) 
S{uf~ - v) 



+ $M1 log 1 _ y *S + S M 1 { log j^z 7 - l°g S ( u l ~ v ) ) + D M\vw ■ 



Simplified equations for Y± 

log|t = \og(l + Y Q )*K Qy (6.17) 

- ^\Og S U y( Ui ,V)+I0g ( , 

i S 2 y{U 2 ' ,v)S 2 y{u\' ,V) 

logF+YL = 2 log 1 * a - log (1 + Yq) *K Q + 2 log(l + Yq) * Kg 1 * s 

1 + *l\w 

i ^ S 2 (u 2 1)+ -v)S 2 (u^--v) | 21og ^(4 1)+ ,.)^(4 2) ',,) ^ 
5 2 (4 2)+ - w) 5 2 (4 1} - - v) S^(u 2 2)+ ,v)S^ ) ~,v) 
-21ogSl* 1 {u 1 ,v)*s + logS 2 (u 1 -v)*s (6.18) 

Simplified TBA equations for Yq 

• Q>4 

• Q = 3 
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(1 + 



log Y 3 = log S(4 2)+ -v)- log Sin®' -v) + log n V , * a • (6-20) 

• Q = 2 

log Y 2 = log -v)~ log 5(4 1} " «) + log TTXXYTTT^ V * . (6-21) 

Hybrid TBA equations for Yq 

Following [27] we introduce a function which combines the terms on the right hand side of 
the hybrid ground state TBA equation 

G Q {v) = -L TBA £ Q + log {1 + Yq,)*{K%Q + 2s nKg-W) (6.22) 
+ 2log(l + Y llvw )*s*K yQ + 2log(l + Y Q 

— l\vw) ~* ° 

- 2 log jj— ^ is * K$ x + log ^*Kq + log (1 - -L) (1 - ^-)*K yQ . 

Then the hybrid TBA equations for Yq read 

clQ ( v m v ) c2Q / (2)- x a2Q , (1)+ x 

logFo( " ) = G ' w - 108 0?U s ^ (, ' 1 - ) +1 ° g ^ 2) (4"-.,) ^(4-,w 

- log Sl2 x ( u i, v ) + 2l °S S( u i ,v) * P .v. k12 x 

-2log S ^ V K Kl2 x ^Q- (6.23) 

s( u y ,v) 

The exact Bethe equations can be found by analytic continuation of e.g. the hybrid 
equations to the string region. In the next appendix, we will consider them for the case 
L = 6. 

Driving terms in the L = 6 case 

The case on which we focus for explicit calculations is L = 6. There, one has that there is 
always exactly one root tm for any M, so that the driving terms take the explicit form 

®M\w = -log5'(r^ / __ 1 - v)S(rJf +1 - v) (6.24) 
®M\vw = -^A/ilog5'(r -v) 
D +x _ = + log S(r i -v) 

D Q= i = +logSl2 x (r ,v) -logS Q (ro - v)S yQ {r^ ,v) + 2logS(r 
Since the twist preserves one supersymmetry, we have [25] 

L TBA = J + 2. (6.25) 
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Driving terms in the L = 10 case 

As another example, we consider a state with L = 10 for which rapidities are outside the 
analyticity strip. One finds that the auxiliary functions Ym\ w and Y- satisfy quantization 
conditions at three distinct (shifted) rapidities r^,r^ , rj^ for any M. As a result, the 
driving terms are now 

3 

M \ W = -^2logS((r < S_ 1 )--v)S((r ( S +1 )--v), (6.26) 
i=i 

3 

E>M\vw = -hn^2logS((r { ^)~ - v) , 

i=l 

3 

D +x _ = J2^S((r?)- -v), 

i=l 

D Q=1 = £ [logS%L (r!P,v) ~ logS Q ((r»)" - t,) S yQ ((r«)~, v)~ 

3 

+2^1o g< s((r«)-, W ) 
i=l 

Furthermore, in this case we have 

Ltba = J- (6.27) 
6.4 Linearized TBA and exact Bethe equations for L = 6 

To find the first perturbative correction to the asymptotic quantization conditions it is 
convenient to expand the TBA system and exact Bethe equations around their asymptotic 
solution. As discussed, this will leave us with three equations (3.11), two of which are 
complex and conjugate to each other, in three real unknowns Su^\ Re(5u^) and Im^u^)- 
These equations are compatible with the quantization of total momentum (3.14). This 
allows one to find a solution for 5u^ by considering one of the two complex exact Bethe 
equations together with (3.14). 
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To this end, we consider the exact Bethe equation for \ that is 
log(-l) = log Y^) = G^) + 21og (1 + Y 1]vw ) * 5 - log '™ 2 m L s^uuvg)) 

^s[(2)( U 3 ' u 2 ) 

- 21o /(^4V g n a _ log5i(r _ _ 41)^^41)) 

c21 C ?/ ( 2 )- V W\ R ps o21 C7/ (1)+ 7/ (1)n | 

^1(2)^3 ' U 2 ) nes ^[(2) l"2 ^2 J 

+ g 021 „(i)a g 021 fJ 2 >+ J 1 )^ w y^( 1 )+^ 1 ' 

^1(2)^3 ' M 2 ) ^1(2)^2 > U 2 jKesy 2 (,^ 2 ) 

+ 21ogResS'*if^(u^,i4 1 " ) ) - log ^ui - 



+ 2 log Res 5 * K^ x (r , u^) - log (V - 




where we used the fact that U2 lies in the overlap of string and mirror regions, and intro- 
duced the short-hand notation 



log 



cy„(2)+ (l)x . c^,( 2 )+ f \ 

W^tt* K ™ 55 / d^og^^KUt,u?). (6.29) 



We now want to expand this and the other TBA equations, by considering 

5(r-M-i-^)'9( r M+i- v ) 
5(rS"-i;) 



Y M |,» = Y^»(l + ^.Wl ^f 1 JeL J . M > 1 (6-30) 



Y Mvw (v) = Y^>)(1 + — r, 

D(r — v) 

Y M \vw(v) = Y^ lv Jv)(l + & M \vw(v)), M > 2 

S(r± — v) 



Y±(v) = YZ(v)(l + & ± (v)) 



S(r{- - v) 



Here Y°(v) are computed out of the asymptotic transfer matrices evaluated at the excep- 
tional rapidities up . These vanish at some root r° = r^uf 1 ) that is not the exact root 
r* dictated by the quantization conditions coming from TBA. The S-matrices on the right 
hand side have poles at these roots, so that the corrections $C are always small on the real 
line. For any Y-function it is convenient to introduce 

A « = rrfe- < 6 - 31 > 

(2) 

Since in many equations terms involving u 2 3 occur, we will have to consider their variation. 

In particular, it is convenient to express them in terms of the difference between u^l and 

(2) 

u 2 3, which we will indicate as 

<5u 2 ,3 = 45 - 45 . ( 6 - 32 ) 
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and for which we know an asymptotic expression (3.4). This quantity should not be 
confused with the corrections Su^l which are the quantities that we are looking for, and 
which cannot be found from asymptotic considerations. In a similar way, we also write 

8r M = r M ~ r° M . (6.33) 

We now proceed expanding the TBA equations. 

Expansion of Y M \ W equations 

%l\w = ^m-i\w [%f-i\w ~ 2ms(r° M _ 2 - v)8r M -2 - 2ms{r° M - v)8r M ) *s 
+ A M+1 \ W {&m+i\w ~ 2m s(r° M - v)5r M - 2m s(r° M+2 - v)5r M +2) * s 

+ Smi + yo&+ ~ (^&- ~ ^) *(rl - v)6r i y j is . (6.34) 

Expansion of Y^ vw equations 

• M = 1 

&i\vw = M\vA\vw * s + (-A_^_ + - 2™s(r? - v){A_ - A+) Sr x ) is 

+ 2m' s(u 2 — v)5u2 — 2iri s(u^ — v)5u 3 — Y 2 ° * s . (6.35) 

• M = 2 

%.\vw = A i\vw^i\vw * s + A 3 \ VW W 3 \ VW -ks-2m A 1][VW s(rl - v) 5r * s 

-Y 3 °*s. (6.36) 

• M > 3 

&M\vw = ^-M-l\vw^M-l\vw * s + Am + i\ vw &m+1\vw * s ~ ^M+l * s ■ (6.37) 

Expansion of Y± equations 



i og y + /y_ 



&+-9L =YZ* K Qy - 2m K 2y (u° 2 ' + , v)5u 2 + 2m K 2y (u° 3 ~ ,v)6u 3 (6.38) 



• iogy+y_ 

W+ + 9L = 2A 1 \ VW [W x \ vw - 2m s(r° - v)8r ) * a 

— 2Ai\ w {&x\ w — 2iri s(tq — v)5tq — 2m s(r 2 — v)5r 2 j * s 
-Y^*K Q + 2Yq * * s + 2m K 2 {u\ - v)5u 2 - 2m K 2 (u^ - v)5u 3 

- Am (K%l(u+,v)5u 2 - K 2 x l v {u 3 ,v)5u 3 ) * s. (6.39) 
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Expansion of the quantization condition for ro 

Since ro appears explicitly in the exact Bethe equation (6.28), it will be necessary to 
consider its quantization condition. The quantization condition for ro should be found by 
continuing the equation for Y_ to —i/g. One can however check that (6.38) is subleading 
in g, so that we can directly work with the equation for log Y+K_ and continue this down 
to —i/g. We have 

log Y+YI = 2 log \ + J?} VW *pv a - log (1 + Y Q ) v K Q (6.40) 

1 + Y l\w 

+ lo S S 2 ( Ul - v) V S + log - log >--±f 

+ log % - « + log "7 ; - - log (i + y^)) , 

where we dropped all the contributions of kernels sub-leading in g. Evaluating this equation 
at ro yields a quantization condition, which can be expanded as follows: 

= 2A 1]vw (&{\ vw - 2m s{r° Q - v)5r ) -k s 

- 2A 1]w {W x \ w - 2m s(rg - v)Sr - 2m s{r° 2 - v)8r 2 ) * a (6.41) 

- y^ v K Q + 2m K 2 {u+ + - r° )5u 2 - 2m K 2 (u 3 - r° )5u 3 - iy^rg) 



+ 8n 



1 + Y l\vw 

i + y° 



2 lo S -i , *P" S ' - lo g *Q *pv ( K n )' + lo § ^(Wl - v) *pi> s' 



+ log S 2 ( U1 - «) v S' + - r °) - ^(n 3 + - r °) + tf(r? - rg) + ^Y^Jrg) 

1 ^"(rg) 1 

2 1 + yi(rg) 2' 



-^i|>g) " o i 1 l ^ + o^i - r° ) 



where the primes denote derivatives with respect to the argument where ro is inserted. 
Expansion of exact Bethe equation for u 2 

From the expansion of the exact Bethe equation we will be able to find the form of 51Z(2), 
as outlined in (3.11). Some care is needed in dealing with the expansion of 

Res^y^,^) _ iog Res5g 2) (4 1)+ ,4 1) ) y^-y^ 
° S5 S2)(4 2)+ ) 4 1) )Resy 2 (4 1)+ ) ° S Res5 s 2 [ 1 (2) (4 2)+ ,4 1) ) ResY 2 (4 1)+ ) 

that according to (3.4) can be written as 

- 2mResK^ {2) (u° 2 + , u° 2 ) 8u 2 - (u° 2 + ) . (6.43) 

The remaining terms can be readily expanded. Since we are interested in the lowest order 
correction to the quantization condition, we can also drop any sub-leading contribution in 
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g, and in particular the terms containing the convolution * . This leaves us with the final 
result 

5Tl {2) = Y%* (k% 2) + 2s * K^ 1 ) + 2A 1]vw (& 1]vw - 2m s(r o ~ - v)Sr ) * S 
- Am (s(u 2 + ,u 2 )5u 2 - s(u° 3 ~ ,u 2 )Su 3 ) * 

+ 4m Ress * Kll w (r^,u° 2 )5r - 2vri/C(rg, u o 2 )5r , (6.44) 
where have introduced the notation JC(u,v) = 2S^5(n, v) with 

S(u,v) = S 1 (u--v)S yl (u-,v) («- -vf ( X i { "\~ X i^\ ) 2 , (6.45) 

\x7(u) - x s (y)J 

in order to conveniently group all driving terms involving ro- 
Cancellation of the most (^-divergent terms 

Finding an explicit expression for the contribution of 5TZ^ to Su 2 and in turn to the 
seven-loop energy is highly non-trivial. The task is much more complicated than in the 
case of the Konishi operator [30, 31] because in (6.44) the correction <5ro appears explicitly, 
together with @{\ vw . To determine these, one would have to solve both linear system 
associated to &m\vw an d to "3^m\vw^ together with the equation yielding the quantization 
condition for 5r$. All these are coupled which makes finding a solution, even numerically, 
a complicated task. 

For the purpose of finding evidence of a non-trivial cancellation of the divergent terms 
in the energy at g 14 , however, a much simpler analysis suffices. 

Let us consider the 0(g 12 ) part of (6.44) and of the linearized TBA equations, and 
expand them in powers of <j). This expansion is expected to involve negative powers, which 
should be the ones that cure the divergences in the energy and that will come multiplying 
the sources of the linear systems. 

For instance, in (6.35) the sources are 

2m s(u 2 — v)5u 2 — 2m s(uj — v)8ii3 — Y 2 * s 
= -2mResY 2 °{u° 2 + ) + 2mResY 2 ° {u°f ) -Y 2 °*s = O(0" 6 ) , (6.46) 

due to the pole of Y 2 at ~ u 2 + + 0((j) e ). This implies that we can expect that 
&x\vw = C^ -6 )- Carrying out a similar analysis for all the remaining TBA equations and 
quantization conditions for auxiliary roots, one concludes that indeed &\\ VV3 = 0((p~ e ) and 
8r M = O(<f> ). 

Turning now to 51Z 2 , we find that up to higher orders in <j> we have 

6K {2) = 7 2 ° * + 2« i^; 2) (,r , «^) fa 3 - (^ + ) + 0(f 6 ) ■ (6.47) 

These three terms are all divergent at 0(0~ 12 ) due to the singularities of Y 2 and K^^y and 
their contribution can be immediately evaluated in terms of asymptotic formulae. Inserting 
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this into (3.11) and using that P = 0((p), one finds indeed that the most divergent part 
of the asymptotic energy at g 14 , which goes like 0((p~ 8 ), is precisely canceled by wrapping 
corrections in (3.5). 
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